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ABSTRACT 

Free electron number densities N e are calculated for cesium plasmas with low electron tem- 
peratures T e {0 . 1 5 - 0. 3 eV) using 5-, 3-, and 2-level atomic models. These Ne values are ob- 
tained for a Maxwellian distribution of free electrons; this distribution was shown to be a good 
approximation for T e j^ 2100° K in optically thick Cs plasmas. The sensitivity of N e to the 
choice of electronic levels, statistical degeneracy and radiative capture coefficients is presented 
for plasmas either optically thin or thick to resonance radiation. The steady state equation for the 
electron number density is discussed in the limits of low and high neutral density Nq s . The ion- 
ization fraction f = N 6 /Nq s was studied as a function of Nq s for the range of T e from 
1740° K (0.15 eV) to 3000° K. Results are presented as a function of T e for N^ s = 3x10^ m"^, 
(typical values for certain thermionic diodes). Both sets of results are compared with the cor- 
responding Saha curves for plasmas in local thermodynamic | equilibrium (LTE). A dis- 

cussion of rules for choosing reasonably accurate but simple atomic models is included. Those 
electronic levels which serve as free electron source terms are identified as functions of Nq s 
for fixed T e . Even a 5-level model gives N e values which are only in fair agreement with 
optically thick results for a 26-level model under certain conditions. The lumped level degen- 
eracy gjyj can be made a function of Nq s and T e to bring the M = 5 results into better agree- 
ment with the more reliable M = 26 results of Norcross and Stone. However, in the absence of 
such a comparison it does not appear that construction of accurate but simple atomic models is a 
straightforward process. 

INTRODUCTION 

There is much interest in the accurate calculation of free electron number density N e in non- 
equilibrium plasmas of low electron temperature T e that exist in Cs plasma diodes. The pur- 
pose of this study is to arrive at rules for constructing simple but satisfactory model atoms. Such 
models must ensure reasonably accurate values of N e for the aforementioned plasmas as well as 
allow simplified calculations of the free electron distribution function. The accuracy of such 
models depends upon the number and arrangement of electronic energy levels as well as inelastic 
cross sections, radiative lifetimes and radiative capture coefficients. For the steady state plas- 
mas of interest, the free electrons are assumed to have a Maxwellian distribution. Solutions of 
the Boltzmann equation (including all elastic and inelastic collision terms) along the lines described 
in Ref. 1 indicate that the distribution function f e (u) is Maxwellian down to T e values of 2100° K 
(for Cs plasmas optically thick to resonance radiation). Although the numerical method described 
in Ref. 1 was convergent only for T e a 2300° K, subsequent calculations down to T e = 2100° K 
have been found to yield essentially Maxwellian f e (u). The high energy tail of f e (u) may be non- 
Maxwellian for T e values <2100° K but the N e values can be relatively insensitive to this tail 
(u > 3. 89 eV) for optically thick plasmas. This insensitivity is present where ground state ion- 
ization becomes a relatively unimportant rate process. 

DETAILS OF THE CALCULATION 

A model atom consisting of M levels is studied. The rate of change of the number of bound 
electrons in any level L is*: 
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N L = -N L (N e X L +^L)+Ne T N K K K _* L + £ N K A K _ L + NeN^K^Ne + /3 L ) (1) 

K?L K>L 


In Eq. (1) Kl-*k (m^ sec -1 ) is the excitation coefficient for the collision induced transition lr*-K 



+ K^ n where K^ n (m 3 sec -1 ) is the ionization coefficient for level L. A k _^l (sec -1 ) is any 
radiative transition probability for de -excitation into the Lth state and is the sum of the 
Aj^k (sec -1 ) for radiative transition to all levels below L. 

For the steady state treated here N L = 0 and the number densities must satisfy the plasma 
normalization condition: 

M 

Ni+X>L = MCS (2) 

1^1 

Equations (1) and (2) together with the condition of charge neutrality (N e = Nj), serve to determine 
N e and the values for a specified f e (u) and initial Cs number density Nq s . Equation (2) can 
be rewritten as: 

M 

f + T\ n l = 1 ^ 

where f is the ionization fraction and is the normalized (to Nq s ) population of the Lth state. 
In the steady state the total ionization rate from bound levels must balance the total capture rate of 
free electrons into all such levels. This relation provides a useful check on the numerical calcu- 
lations. Its application indicated that accurate numerical solutions for f were not possible at low 
values of N e withput tight convergence criteria. 

The condition, N e = 0 can be written as: 

M M 

[(Kf P)*f 2 + 0 L f] = ^ Kj? n Ni (3) 

where (K ca P)* (m 3 sec -1 ) is an effective two-body capture coefficient = N^gK^^. Thus the 
equation for the ionization fraction is essentially quadratic; however, the normalized populations 
Nl are strong functions of Nq s and optical thickness for fixed T e . The ionization coefficient 
K *£ n — (Q ion (v e )v e ) is of course a function of T e alone. 

MODEL ATOM PROPERTIES 

The optimization of the model atom depends upon the choices of cross sections for electron im- 
pact and coefficients and frequencies for radiative processes involving the M*h level. This level 
must be adjusted so as to simulate the presence of many missing excited levels 1 . The assign- 
ments of ionization potential Ej^ statistical degeneracy gj^ and Jlyi and /3j^ values thus be- 
come critical for accurate calculations of f values. The Nj^ for this lumped level is a com- 
plicated function of the level structure for variable Nq s (see Eq. (1)). The schematic diagram 
for the three models is shown in Fig. 1 with values of E L , Al-^k SL* 

The early form of the Gryzinski cross sections has been used to compute ionization coefficients 
and excitation coefficients for optically allowed transitions only. 2 The Gryzinski exchange formu- 
lation was used to compute excitation coefficients for optically forbidden transitions. 3 Although 



the exchange cross sections have relatively small values at maximum, their slopes are steep just 
above threshold. Since the excitation and ionization thresholds are usually »kT e , the behavior of 
the cross section in this region determines the value of the excitation coefficient. The lumped 
level was treated as an optically allowed state for all excitation (and de -excitation) collisions. 

Most allowed and exchange cross sections (monoenergetic) give reasonable agreement with ex- 
periment over a range of electron energies from threshold to maximum 4 . The experimental value 
of the most important excitation cross section Q®* 2 ^ as a slope roughly 3 times the Gryzinski 
value at threshold 4 . 

The radiative capture coefficients /?l were taken from Ref. 5 where they were calculated 
using an adjusted quantum defect method. These calculations agree well with known oscillator 
strengths and recombination cross sections. 

The radiative transition probabilities were taken from Ref. 6; the agreement of those calcula- 
tions with cited experimental values is within 50% for all important A LK values (see Ref. 7). 

The most critical Alk value is for the 2->T resonance transition; this line is strongly absorbed 
for all experimental plasmas of interest 8 . For M = 3 the A32 value was set equal to the sum of 
A52, A42 3 and A32. The value of A21 was used for all three optically thin cases. 

RESULTS - OPTICALLY THIN PLASMAS 

Since the only available M - 26 results (optically thin) are several points for T e = 3000° K, 
the M = 2 and 3 results will be compared with 5-level results for those plasmas. Since ground 
state ionization is a significant process at low Nq s values, however, a non-Maxwellian tail in 
f e (u) could cause large changes in calculated values of f in this region. 

Low N° Cs limit. - In the limit of low Nq s values three-body recombination is negligible and 

the free electron number density is determined by (from Eq. (3)) 

M M 

E *L fs E K L nN L (4) 

L=1 L=1 

hi the optically thin case, radiative capture into excited states is the important capture process 
which balances ground state ionization. Then f is given approximately by: 


f « Kj° n N. 



(4a) 


which simplifies for M = 2. The product Kj° n N^= Kj° n (for allM), however, since the ground 
state population is greater than 0.95 Nq s for optically thin plasmas in the low Nq s regime. The 
value of f is then simply determined by the ratio of the ground state ionization coefficient to the 
sum of the largest radiative capture coefficients. For cesium the /3g (5D level) value dominates 

the Pl term; it is three times the maximum / 3 value for other excited levels, /3(6P), and 
200 times the (6S) value 5 . Thus, in choosing simple model atoms, it is important to assign 
a large value of to simulate the missing excited levels at low The /3 5 assigned for 

M = 5 is the sum of |3 values for the six excited states: 7P, 6D, 8S, 4F, 8P, 7D. 

The results of f versus Nq s for 2-, 3-, and 5-level models are shown in Fig. 2 for an elec- 
tron temperature of 3000° K (0.26 eV). The f values for M = 2 (6S and 6P levels) with /? 2 = 
/3(6P) + /3(5D) agree with the M = 5 results within 50% at low values of Nq s (< 3xl0^ 5 m“3). 

This agreement can be improved by adding the sum (/3(7S) + /3g) to the /3 2 value. For a given 
excitation cross section the de-excitation coefficient is inversely proportional to 

the upper state degeneracy gj£. The L^ 1 level 3-body capture coefficient is directly propor- 
tional to the gL value. The change of degeneracy for the M^ 1 state (M - 2,3) does not affect 
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the f value at low Nq s because the value of is unchanged. However, the degeneracy does 
affect the maximum ionization fraction calculated for the M = 2 and 3 models. This latter be- 
havior can be described in terms of the free electron source terms although there is no simple 
limiting expression. The degree of departure of f from LTE conditions is evident by direct 
comparison with the Saha curve in Fig. 2. 

High Nq s limit. - In the limit of high values 3-body capture becomes the dominant loss 

mechanism for free electrons. Radiative capture is negligible and the value of f can be approxi- 
mated by: 



The excited levels have the largest ionization coefficients since the ionization cross section Ql 
is a strong inverse function of ionization potential El. These levels are most nearly in equilib- 
rium with the free electrons and f approaches the Saha result. For the simplest model atom, 

M = 2, one must choose values of g M and E]y[ so as to preserve the ’'correct” ratio of ioniza- 
tion rate to 3-body capture coefficient. In this study the value of E 2 was fixed at the 6P value, 
2.43 eV, and the value of g 2 was varied for best agreement with M = 5 results. However, it 
will be shown that the M^h state for M = 5 is not so simply adjusted to match the M = 26 results. 

For a given ratio of ionization and capture coefficients for the lumped level N 2 must vary with 
NCS so as to account for the populations of the missing excited levels near the Saha limit. For 
M = 2, N 2 increases so rapidly (cc Nq s ^“ that K 2 011 exceeds (K ^ a Pf*f 2 at intermediate 
values of Nq s . The f is relatively small and only a gradual function of Nq s so the 3-body rate 
does not become important for the optically thin case until high Nq s values. The maxima in the 
f curves of Fig. 2 occur approximately where the collisional lifetime = ^Ne)*"^ e( l ua l s the radi- 
ative lifetime (A 2 i)"^ for the first excited state (see Ref. 9). 

The M = 3 and 5 cases are similar in that net ionization from excited states is a significant 
process for Nq s > 10 22 m -2 . The M = 2 model is unique because the ground state always serves 
as the free electron source term (i.e. exhibits net ionization). Thus in this case the lumped level 
necessarily becomes a net capture r since N e = 0. The large ( 3 2 value (60 times /3j_) ensures 
that this level remain a net capturer at low values of Nq s . The sources of ionization for the 
M = 3 and 5 models can be determined from the net ionization rates. The rate for the L& state 
can be written from Eq. (3) as: 

R net _ K ipn _ j(KffP)*f 2 + (6) 

Only R^ et values which are positive will be considered, i.e., states which are free electron 
source terms. The variation of such source terms (normalized to the sum of positive terms) 
r L is shown for T e = 3000° K, optically thin, in Fig. 3 for M = 5 and 3. The L = 4 level 
(7S) makes the largest contribution to the free electron population at T e = 3000° K for M = 5, 

Nq s > 2xl0 21 m“ 2 . This behavior differs slightly from the T e = 0. 2 eV results where the L = 3 
state (5D) has the largest r£, value at high Nq s . 

The results for M = 5, g 5 = 50, Eg = 0. 6 eV of Fig. 2 are in good agreement (within a factor 
of two) with Ng/ifNgJgaka results for the 2 6 -level model at T e = 3000° K 8 (see Table I). No 
optically thin results for M - 26 were available for comparison at lower T e values. The sim- 
pler models should give good agreement with the M = 26 case at low N^ s (<10 19 m -2 } since the 
lumped /3-jyj- values have been chosen so as to produce that behavior 

RESULTS - OPTICALLY THICK PLASMAS 

Low N^ s limit. - At low values of Nq s ( <10*9 m -2 in Fig. 4(a)) f is determined by the bal- 
ance between ionization from excited states and radiative capture. Just as in the optically thin 
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case for high N^g, the excited states are free electron source terms, but at low Nq s . | The 
dominant is R 2 for Nq s ^ m“ 2 in the 3- and 5-level models. Only the ground state 

has a positive R^ et value for M - 2 through the range of Nq s values. The approximate rela- 
tion for the ionization fraction becomes: 



For M = 5 the above ratio should give reasonably accurate f values when just the L = 2 and 3 
ionization terms and the /3g value are included. However, at intermediate values of Nq s 3-body 
capture is comparable with the radiative process and the lumped level must represent the missing 
excited states. The excited states are free electron source terms at low Nq s because their 
populations (especially N 2 ) are enhanced in the absence of radiative de-excitation. The ioniza- 
tion rates become large even when populations Nl, L > 1 remain & 10" 2 ■-> 10“ 2 of because 
the ionization coefficients K^ n » K^ on . The mono energetic cross sections Q^ n , L > 1, have 
lower energy thresholds and much higher slopes than Q| on in the threshold region. The excitation 
coefficients for collision-induced transitions between excited levels are correspondingly larger 
than the ionization coefficients and have low thresholds (some <1 eV). 

The f values for the optically thick case are plotted versus Nq s for T e - 3000° K in 
Fig. 4(a). The lumped /Jjyj value should guarantee accuracy of the simpler models at low Nq s . 
This is true if the populations N 2 have nearly the same values as for M = 26. Decreasing /%[ 
for M = 2 gives better agreement with M = 26 at Nq s = 10 2 ^ m” 3 but the asymptote is incor- 
rect at 10^® m -2 . For M = 3, the f value is increased at all N^ s by making the plasma thick 
to all line radiation and lowering g^ from 50 to 10. There is no peak in f versus Nq s for 
M = 2, g 2 = 50 since the L = 2 state has a large 3 -body capture rate at intermediate Nq s due 
to its large g^ value and relatively high N e . Lowering the degeneracy to g 2 = 6 (not shown) 
does not change the curve shape since 3 -body capture continues to dominate excited state ioniza- 
tion. Decreasing the g^ from 50 (curve 4) to 10 (not. shown) for M = 3 does produce a maxi- 
mum in f because 3->2 de-excitation is increased enough to raise the n}j value significantly. 
This L = 2 state is the main free electron source term for M = 3. 

The f values for the optically thick case, T e = 2321° K, are plotted versus Nq s in Fig. 4(b). 
Just as for T e = 3000° K the f is determined by ionization from the ground level for M = 2 
and from excited states for M = 3 and 5. These ionization processes balance radiative capture 
into all levels and the curves remain flat where the populations of the ionization sources remain 
constant, just as for 3000° K. The features of the curves are identical to the T e = 3000° K re- 
sults throughout the range of Nq s values. The peak for M = 5 is again caused by ionization 
from the higher excited states L = 3 and 4. Four points of M = 26 results are also shown in 
Fig. 4(b). Comparison with the latter results indicate that the simpler atomic models overesti- 
mate departures from LTE at values of Nq s below 3xl0 22 m“ 2 . It should be noted, however, 
that the N e values from M = 26 are considerably lower than Saha values for Nq s as 10 2 ^ m” 2 . 
The ionization fractions are approximately 10 -2 of the 3000° K values because excited state 
populations and ionization coefficients are considerably lower. The limiting (low N^, g ) value of 
f for M = 3 is only 10 ” 2 of the f computed for a 3 -level model in Ref. 9. For M = 2 the value 
of f is about 1/15 of that for a 2 -level model in Ref . 9. These large differences are due to the 
present use of more reliable excitation cross sections and and values. 

The variation of the free electron source terms for M = 5 at T e = 2321° K (0. 2 eV) and 
3000° K are shown in Figs. 5(a) and 5(b), respectively. Qualitatively, the dependences of Rp, 
are similar but the region of rapidly changing contribution is shifted by a factor of 10 2 to higher 
Nqs values at the lower value of T e . The M = 26 results of Ref . 7 indicate that the maximum f 
is nearly 5 times the M = 5 value for N^ g = 10 2 ^ m” 3 , T e = 2321° K. Figure 5(b) shows that 
the main contributor of free electrons is the L = 3 state for Nq s > 3xl0 20 m” 2 and the ground 
and first excited states for Nq s values ^3xl0 20 m” 2 . If additional excited states are to become 
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new ionization source terms, their relative populations must fall off slowly enough with binding 
energy that their higher ionization coefficients can raise f , The sudden drop in the L = 3 contri- 
bution in the region Nq s s lO 2 -^ m -3 is a result of the large /3g value since radiative capture 
becomes the dominant electron loss process in that region (N e sS SxlO* 8 m -3 ). 

The agreement between the M = 5 and 26 results improve at high Nq s (= 10 22 m -3 ) for 
T e = 3000° K as can be seen in Fig. 4(a). The f for M = 5 is about 90% of the f value for 
M = 26. Also, the disparity near peak (N^ s = 10 20 m -3 ) is explainable in terms of the value 
just as the T e = 2321° K results. Figure 5(b) shows that the L = 4 level is the chief contributor 
of free electrons for Nq s > 3xlo 19 m -3 . This occurs because the "ladder" ionization mechanism 
operates if the excited state populations are large enough. For optically thick plasmas this is the 
case. It is plausible, therefore, that including additional excited states would raise the value of f 
in this T e region if these states serve as large source terms with R L values = RV. 

High Nq s limit. - At high values of Nq s the excited levels are still the free electron source 

terms. Again the balance which determines f is between excited state ionization and three-body 
capture into these same levels. The f values, optically thick, are uniformly larger at lower 
N^g values than in the optically thin case. This is because the important Nl approach their 
Saha values at lower Nq s (s 3x 1 m -3 ). The value of f can be satisfactorily calculated from 

Eq. (5) including terms, L-2 to M. 

Results for T e = 1740° K (0. 15 eV). - Curves of f versus Nq s were also calculated at 
T e = 1740° K, There is no structure in these curves for N^ g < 10 23 m -3 , optically thin and 
Nq s < 10 28 m -3 , optically thick. One point was compared with the M = 26 results at Nq s = 

10 24 m -3 , optically thick. The M - 26 model gives an f which is 3 times the value for M = 5, 
g5 = 50. Since the chief free electron source term is L = 5 for M = 5 the difference between the 
M = 5 and 26 results becomes very sensitive to the lumped state. This level must now simulate 
the source terms for N e . 

The simpler models do not agree well with each other at T e = 0. 15 eV, optically thick. This 
behavior is shown in Fig. 6 where f is plotted as a function of T e at Nq s = 3xl0 22 m -3 for 
M = 2, 3, and 5. The largest disparity occurs in the optically thick case where the f for M = 2 
is only 2% of the value calculated for M = 5, This result is consistent with the fact that the 
former model does not allow for excited state source terms. The M = 3 model gives an f value 
midway between the M = 2 and 5 results. The values of f, optically thin, are only = 10 -3 of the 
optically thick values; this is because the ground state is the sole ionization term in the former 
case. Thus at this low T e value the ionization mechanism changes from simple ground state ion- 
ization to a "stepwise" ladder process (for Nq s = 10 18 - > 10 25 m“ 3 ) as the optical thickness 
changes. As the electron temperature is increased this disparity between optically thin and thick 
results decreases markedly. For fixed values of /3|^ and E^, agreement with M = 26 results 
for simpler models can be improved by using variable gj^ = F(Nq s , T e ) or making the plasma 
thick to all radiation. However, any scaling laws for construction of atomic models remain, to be 
investigated in the absence of more detailed results for comparison. 


CONCLUDING REMARKS 

For certain values of electron temperature and neutral density 2- and 3-level models can be ad- 
justed to give N e values which agree with results for 5-levels. However, for the experimentally 
interesting case of optically thick Cs plasmas, results for a 26 level model 7 indicate the "best" 

M = 5 model can seriously underestimate the ionization fraction f . This is especially true for 
values of electron temperature <3000° K. This disparity occurs because the sources of free elec- 
trons are highly excited levels as a result of a "ladder" ionization mechanism. It does not appear 
likely that convenient adjustment of simpler models will improve agreement with the M = 26 re- 
sults. For proper choices of E]^, for the lumped level a variable gMC^Cs) wi ^ Pro- 

vide satisfactory agreement for fixed T e . However, it appears that increasing the number of 



levels is a more desirable approach than this to ensure calculation of accurate Ne values. Fur- 
ther results for M = 26 are needed to evaluate the simpler model atoms for the low temperature 
(1740° - 2321° K) optically thin plasmas. 


TABLE I. - (T e = 3000° K; OPTICALLY THIN) 


M 

Ne/ (N e )saha 

5 

2.4X10- 2 

6.8x10-2 

0.35 

0.75 

26* (Ref. 7) 

=1. 3x10 ”2 

=6.5x10-2 

= .32 

= .75 


♦Points estimated from plots without grid. 
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Figure 1. - Schematic diagram of cesium model atoms. 



Figure 3. - Plots of normalized net ionization 
rates as a function of Cs neutral density in 
plasmas optically thin to resonance radiation 
for 3- and 5-level atomic models. T e - 3000° K. 



Figure 2. - Plots of ionization fraction as a func- 
tion of Cs neutral density in plasmas optically 
thin to resonance radiation for 2-, 3-, and 
5-level atomic models. T e = 3000° K. 


(a) T e = 3000° K. Two points for a 26-level 
model are shown for comparison. 

Figure 4. - Plots of ionization fraction as a func- 
tion of Cs neutral density in plasmas optically 
thick to resonance radiation for 2- 3-, and 
5-level atomic models. 
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(b) T e = 2321° K. Four points of 26-level results 
are included for comparison. 

Figure 4. - Concluded. 



(a) T e = 3000° K. 


Figure 5. - Plots of normalized net ionization 
rates as a function of Cs neutral density in 
plasmas optically thick to resonance radiation 
for the 5- level atomic model. 



Cs number density, Nq s , m‘ 


(b) T e * 2321° K (0. 2 eV). 
Figure 5. - Concluded. 
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Figure 6. - Plots of ionization fraction as 
a function of electron temperature in 
plasmas optically thin and thick to reso- 
nance radiation for M = 2-, 3-, and 
5-level atomic models. 



